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Abstract 


A  computer  is  used  to  simulate  a  co.lisionless  plasma  in  an 
axisyanetric  high  specific  impulse  plasma  accelerator  with  a  constant  . 
applied  electric  field  and  with  no  applied  magnetic  field.  The 

collisionless  plasma  model  is  shown  to  be  applicable  to  plasma  acceler- 

k 

ators  that  operate  with  specific  impulses  in  excess  of  10  seconds  and 

15  ^ 

with  particle  number  densities  of  the  order  of  10  '  particles/cmJ  or 

les^.  In  particular,  the  model  is  used  to  simulate  an  accelerator  that 

has  recently  produced  specific  impulses  of  the  order  of  10^  seconds 

15  ^ 

with  particle  number  densities  of  the  order  of  10  particles/cm  . 

—  The  computer  simulation  is  accomplished  by  representing  the  many 
ions  and  electrons  in  the  plasma  by  several  thousand  representative 
particles.  These  several  thousand  particles  are  assumed  to  constitute 
a  good  statistical  sample  of  the  actual  ions  and  electrons,  and  the 
normalized  distribution  function  of  the  plasma  is  approximated  by  the 
normalized  distribution  function  of  the  representative  particles.  *  The 


motion  of  these  representative  particles  is  follpJBBd-in  tinurTand  the 
electromagnetic  fields  induced  by  these  particle  motions  are  calculated 
by  integrating  the  charge  and  current  densities  obtained  from  the 
particle  distributions. 

S  A  computer  program  has  been  developed  to  display  simultaneously 
on  a  cathode  ray  tube  the  electric  potential,  the  magnetic  field,  and 
the  motion  of  the  representative  particles  as  they  move  through  the 
accelerator  under  the  influence  of  the  fields.  ^>Xbese  displays,  which 
have  been  recorded  on  motion  picture  film,  illustrate  the  coupling 
between  the  particle  motions  and  ths  fields.  In  particular,  the  dis- 


iii 


plays  Illustrate  how  the  plasma  motion  modifies  the  applied  electric 
potential,  how  the  currents  induce  large  magnetic  fields,  how  the  ions 
receive  their  energy  directly  from  the  electric  field,  and  how  the 
ions  are  deflected  toward  the  accelerator  exhaust  by  the  induced 
magnetic  field. 

The  results  of  the  computer  simulation  illustrate  that  two  fre¬ 
quently  used  concepts  for  a  nearly  neutral  highly  conducting  plasma 
should  not  be  applied  to  this  problem.  These  concepts  are  (l)  that  the 
PjK  force  cancels  the  p@E  force  for  a  neutral  plasma,  and  (2)  that  ^ 
plasmas  with  infinite  conductivity  are  frozen  to  magnetic  field  lines. 
Although  the  plasma  is  nearly  neutral,  the  p^E  force  is  not  strictly 
cancelled  by  the  p  E  force  as  in  a  highly  collisional  plasma  because 
the  electron  and  ion  gas«:s  can  easily  slip  through  one  another.  For 
this  case,  the  p^E  force  on  the  electrons  is  balanced  almost  entirely 

by  the  v  xB  force  rather  than  by  collisions  with  the  ions,  and  the 
© 

p  E  force  is  the  dominant  force  for  adding  energy  to  the  ions  just  as 
in  the  electrostatic  ion  gun.  Although  the  electrical  conductivity 
in  the  simulation  is  effectively  infinite  because  there  are  no  particle- 
particle  collisions  included  in  the  simulation,  the  fluid  is  not  frozen 
to  the  magnetic  field  lines  because  oi  the  large  ion  inertia  forces. 

Although  the  ions  are  accelerated  by  an  electrostatic  field  just 
as  in  an  ion  gun,  the  current  is  not  space -charge  limited  as  in  an  ion 

i 

gun  because  electrons  are  present  to  partially  neutralize  the  plasma. 

The  electrons  are  restricted  from  flowing  directly  to  the  anode  by  the 
self-induced  magnetic  field. 

Since  the  ions  receive  their  energy  from  an  elr>ctrostatic  field, 
they  do  not  obtain  the  "anomalous  velocities"  (i.e.,  velocities  greater 
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than  they  would  receive  by  falling  through  the  potential  applied  across 
the  electrodes)  that  are  reported  for  the  MPD  arcjets. 

Performance  predictions  based  on  the  computer  simulation  are 
given  for  the  high  specific  impulse  accelerators  for  applied  Voltages 
from  10  to  10  volts. 

The  results  from  the  simulation  suggest  that  the  high  specific 
impulse  plasma  accelerators  should  operate  most  efficiently  with  no 
applied  magnetic  field  and  with  a  mass  flow  rate  proportional  to  the 
current.  The  specific  impul.se  should  increase  in  proportion  to  the 
current,  and  the  voltage  and  thrust  should  increase  with  the  square  of 


the  current 
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I.  INTRODUCTION 


Electric  plasma  accelerators  are  attractive  for  deep  space 
propulsion  because  they  offer  a  required  combination  of  high  specific 
impulse  and  high  thrust  density.  The  high  specific  impulse  is  re¬ 
quired  for  efficient  use  of  the  propellant,  and  the  high  thrust  density 
is  required  to  keep  mission  times  within  reasonable  bounds. 

A  graph  produced  by  Clark  and  Jahn1  that  compares  the  specific 
impulses  and  thrust  densities  of  electric  accelerators  with  specific 
impulses  and  thrust  densities  of  other  propulsion  devices  is  shown  on 

figure  1.  This  graph  has  been  modified  to  show  the  new  state  of  the 

2 

art  for  pulsed  plasma  guns  as  given  by  Cheng  .  The  graph  illustrates 
that  other  propulsion  devices  can  produce  either  high  specific  impulse 
or  high  thrust  density,  but  not  both. 

The  chemical  rocket  can  produce  high  thrust  density  but  not 
high  specific  impulse.  It  converts  chemical  energy  to  thermal  energy 
and  then  converts  themal  energy  to  directed  energy  with  a  nozzle. 

With  this  arrangement,  the  directed  energy  per  unit  mass  of  propellant 
cannot  exceed  the  chemical  energy  per  unit  mass  available  for  conversion 
The  chemical  rocket  can  achieve  a  high  thrust  density  by  using  a  large 
mass  flow  of  propellant,  but  the  velocity  and  the  specific  impulse 
are  limited  as  shown  on  figure  1.  On  the  other  hand,  the  ion  engine 
can  produce  high  specific  impulse  but  cannot  produce  high  thrust  density 
The  flow  of  ions  that  are  accelerated  electrostatically  by  the  high 
voltage  acceleration  grid  is  limited  by  the  space  charge  of  the  ions 
between  the  ion  source  and  the  acceleration  grid.  Therefore,  even 
though  the  ion  engine  can  accelerate  each  ion  to  high  velocity,  it  can 
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not  produce  high  thrust  density.  The  nuclear  rocket,  which  converts 
nuclear  energy  to  thermal  energy  and  then  thermal  energy  to  directed 
energy  with  a  nozzle,  can  produce  high  thrust  density  but  is  limited 
in  specific  impulse  by  the  properties  of  the  material  containing  the 
heated  propellant .  Although  a  large  amount  of  energy  per  unit  of  mass 
can  be  converted  from  nuclear  energy  to  thermal  energy,  the  propellant 
must  be  contained  by  the  physical  wall;  so  the  thermal  energy  is  limited, 
and  consequently  the  specific  impulse  is  limited  to  the  value  shown  on 
figure  1.  In  a  resistojet,  an  electrically  conducting  material  is 
heated  electrically  and  the  propellant  is  in  turn  heated  by  the  con¬ 
ducting  material.  The  propellant  temperature  is  always  less  than  the 
temperature  of  the  conducting  wall,  and  consequently  the  specific 
Impulse  is  limited  to  the  values  shown  on  figure  1. 

The  remaining  devices  shown  on  figure  1  are  electric  plasma 
accelerators.  The  conventional  arcjet  uses  an  electric  arc  to  heat 
the  gas  and  then  converts  the  thermal  energy  to  directed  energy  with  a 
nozzle.  These  arcjets  are  limited  in  specific  impulse  by  the  proper¬ 
ties  of  the  material  containing  the  heated  propellant  exactly  the  same 
as  the  nuclear  rocket.  The  cross-field,  HPD,  and  Cheng's  accelerators 
can  achieve  much  higher  specific  impulses  because  the  electrical  energy 
is  not  all  converted  to  thermal  energy  —  rather  much  of  the  energy  is 
converted  straightway  into  directed  energy,  and,  consequently,  no 
nozzle  is  required  to  convert  thermal  energy  into  directed  energy.  The 
specific  impulse  achievable  is  then  limited  by  the  amount  of  electrical 
energy  that  can  be  converted  straightway  into  directed  energy  without 
impinging  excessive  energy  onto  the  material  components  of  the  acceler- 
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ator.  Some  of  the  propellant  impinges  on  the  accelerator  wall  and 
the  electrical  current  carried  by  the  electrodes  imparts  heat  to  the 
electrodes;  so  the  amount  of  energy  that  can  be  imparted  to  the  gas 
will  ultimately  be  limited  even  in  the  devices  that  convert  a  high 
fraction  of  tue  total  electrical  energy  directly  into  directed  energy. 
This  limit  has  not  yet  been  accurately  determined.  Nevertheless,  it 
has  been  shown  experimentally  that  the  combination  of  specific  impulse 
and  thrust  density  already  obtainable  with  the  electric  plasma  ac- 

i 

celerators  is  much  greater  than  that  obtainable  by  the  other  propulsion 
devices . 

Plasma  accelerators  also  have  an  application  as  a  source  for  high 
energy  plasma  for  controlled  nuclear  fusion  devices.  All  the  plasma 
containment  devices  designed  thus  far  for  continuous  operation  lose 
plasma  at  a  fairly  rapid  rate,  and  in  order  to  replace  the  lost  plasma 
a  source  of  high  energy  plasma  is  required.  Plasma  accelerators  have 
the  capability  to  produce  a  plasma  with  the  energy  required  for  fusion 
at  mass  flow  rates  orders  of  magnitude  larger  than  the  ion  guns  presently 
being  used. 

The  development  of  the  electric  plasma  accelerator  has  been 
primarily  achieved  by  experimental  trial  and  error  techniques  because 
a  good,  detailed  theoretical  analysis  was  not  available.  At  the  end 
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of  1963,  DucatiJ  found  experimentally,  to  his  surprise,  that  the 
specific  impulse  within  an  arejet,  such  as  the  one  shown  in  figure  2, 
could  be  increased  from  the  region  on  figure  1  denoted  as  the  conven¬ 
tional  arejet  to  the  region  denoted  as  the  MPD  arejet  by  simply  re¬ 
ducing  the  pressure  in  the  accelerator  from  near  atmospheric  pressure 
to  pressure  under  0.1  atmosphere,  and  by  increasing  the  current  from 
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hundreds  of  amperes  to  thousands  of  amperes.  In  1968,  Cheng  obtained 

the  further  increased  performance  shown  on  figure  1  from  an  experiment 

in  which  he  maintained  a  low  pressure  and  further  increased  the  current 

from  thousands  of  amperes  to  the  order  of  one-hundred  thousand  amperes 

4 

within  a  coaxial  configuration  as  shown  on  figure  3-  Clark  and  Jahn 
are  experimenting  with  an  electric  accelerator  with  the  configuration 
shown  on  figure  4  at  currents  in  excess  of  100, (XX)  amperes,  but  they 
have  not  yet  measured  the  specific  impulse  and  thrust  density  within 
their  accelerator.  Although  Cheng's  and  Clark  and  Jahn's  accelerators 
are  pulsed,  the  characteristics  of  the  accelerators  are  nearly  steady 

4 

over  most  of  the  pulse,  and  Clark  and  Jahn  point  out  that  the  perfor¬ 
mance  is  a  good  simulation  of  steady  state  accelerators. 

Disguised  by  the  seemingly  continuous,  smooth  transition  of  the 
experimental  results  from  the  conventional  arcjet  regime  to  the  MPD 
arcjet  regime  and  onto  the  regime  of  Cheng's  data  is  the  fact  that  the 
relative  importance  of  the  various  accelerating  mechanisms  differ 
appreciably  in  the  various  regimes.  Theoretical  models  have  recently 
been  developed  to  describe  the  acceleration  mechanisms  for  the  conven¬ 
tional  arcjet  and  for  the  MPD  arcjet.  The  purposes  of  this  dissertation 
are  (l)  to  present  a  theoretical  investigation  of  the  acceleration 
mechanisms  for  an  electric  accelerator  in  the  high  performance  regime 
of  Cheng's  data,  (2)  to  use  the  results  of  this  investigation  to  argue 
that  the  models  recently  developed  to  approximate  the  acceleration 
mechanisms  for  the  MPD  arcjet  regime  cannot  be  extended  into  the  high 
performance  regime,  (3)  to  suggest  a  more  appropriate  theoretical 
model  to  be  used  for  this  regime,  and  (4)  to  present  results  of  a 
computer  simulation  based  on  this  model. 
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II.  ACCELERATING  MECHANISMS  WITHIN  ELECTRIC  PLASMA  ACCELERATORS 


The  shift  in  importance  of  the  different  ace ole rating  mechanisms 
can  be  related  to  the  importance  of  particle  collisions  within  the 
accelerator.  When  the  particle  free  paths  are  sufficiently  small,  as 
in  the  conventional  arc jets  that  operate  at  near  atmospheric  pressure, 
the  electrical  energy  given  to  the  charged  particles  is  quickly  random¬ 
ized  by  collisions.  No  appreciable  directed  energy  is  imparted  to  the 
gas  by  the  electromagnetic  field  —  instead  the  thermal  energy  is  con¬ 
verted  to  directed  energy  with  a  nozzle. 

As  the  charged  particle  path  lengths  between  encounters  increase, 
the  charged  particles  are  deflected  between  collisions  by  the  magnetic 
field  and  receive  a  component  of  velocity  in  the  JxB  direction.  If 
collisions  are  sufficiently  frequent,  most  of  the  electrical  energy  is 
given  to  the  electrons,  and  the  electrons  transfer  their  directed 
energy  to  the  ions  by  collisions.  However,  if  collisions  are  sufficiently 
rare,  the  electrons  cannot  transfer  their  energy  to  the  ions  by  close 
collisions  and  the  ions  receive  most  of  their  energy  directly  from  the 
collective  electromagnetic  fields. 

To  determine  whether  or  net  the  electrons  can  transfer  their 
directed  energy  to  the  ions  by  close  collisions,  one  must  compare  the 
electron  path  lengths  required  to  transfer  energy  to  the  ions  with  the 
dimensions  of  the'  accelerators,  (This  path  length  is  much  longer  than 
a  "free  path"  because  the  energy  exchange  in  an  elastic  collision  be¬ 
tween  a  light  particle  and  a  heavy  particle  is  small  compared  to  the 
kinetic  energy  of  ths  particles).  Accelerator  dimensions  are  usually 
of  the  order  of  centimeters.  The  MPD  arejet,  Cheng's  accelerator,  and 
Clark  and  Jahn's  accelerator  operate  with  ion  number  densities  of  the 
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order  of  10  ions/cm  throughout  most  of  the  accelerator.  For  this 

density,  the  minimum  electron  path  length  required  to  transfer  to 

hydrogen  ions  energies  corresponding  to  specific  impulses  from  10^  to 

5 

10  seconds  have  been  estimated  and  are  shown  in  figure  5«  These  mini¬ 
mum  path  lengths  were  calculated  using  the  coulomb  cross  section  and 
assuming  that  the  average  velocity  of  electrons  is  equal  to  the  average 
ion  exit  velocity.  Electron  velocities  are  most  likely  higher  than  this 
minimum  value;  therefore,  the  corresponding  path  lengths  are  most  likely 
longer  than  the  minimum  values  shown  in  figure  5«  The  path  length 
increases  rapidly  with  increasing  energies  because  the  coulomb  cross 
section  decreases  rapidly  with  increasing  energy.  Figure  5  shows  that 
for  an  electron  to  transfer  10  ov  to  the  ions  a  path  length  of  the  order 

6  k 

of  1  cm  is  required,  whereas  a  10  cm  path  is  required  to  transfer  10 

6V. 

The  range  of  specific  Impulses  depleted  by  Clark  and  Jahn  for  the 
MPD  arejet  is  also  shown  on  figure  5?  and  it  is  seen  that  fcne  corresponding 
path  lengths  are  of  the  order  of  or  smaller  than  the  dimensions  of  a 
typical  accelerator.  Therefore,  in  such  a  device  the  electrons  can 
transfer  to  the  ions  by  collisions  much  of  the  energy  directed  in  the 
JxB  direction.  When,  in  addition,  the  number  of  electrons  that  pass 
through  the  accelerator  is  large  compared  *0  the  number  of  ions  that 
pass  through  the  accelerator  (electron  current  much  greater  than  ion 
current)  each  ion  can  be  given  energy  from  many  electrons.  Thus,  the 
ion  can  achieve  an  "anomalous  velocity",  i.e.  a  velocity  greater  than 
it  would  receive  by  falling  freely  through  the  available  potential 
field. 
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On  the  other  hand,  in  Cheng's  experiment,  the  specific  impulse 
is  of  the  order  of  105  secs  and  the  electron  path  length  shown  on 
figure  5  for  this  specific  Impulse  is  orders  of  magnitude  larger  than 
the  largest  dimension  of  his  accelerator.  Therefore,  either  (l)  many 
more  electrons  than  ions  flow  so  even  though  each  electron  gives  only 
a  small  increment  of  energy  to  an  ion  sufficient  electrons  are  present 
to  give  the  relatively  few  ions  their  final  energy;  or  (2)  the  ions 
receive  most  of  their  energy  directly  from  the  electomagnetic  fields. 

If  the  former  were  correct,  most  of  the  electrons  would  go  to  the  anode 
because,  for  steady  operation  of  an  electrically  floating  accelerator, 
the  number  of  electrons  leaving  the  accelerator  must  equal  the  number 
of  ions  leaving  the  accelerator.  Since  most  of  the  energy  of  the 
electron  could  not  be  transferred  to  the  ions,  most  of  the  energy  would 
be  transferred  to  the  anode  resulting  in  efficiencies  orders  of  magni- 
tude  less  than  unity.  Although  efficiencies  have  not  been  accurately 
measured  in  Cheng's  accelerator,  order  of  magnitude  measurements  shew 
that  they  are  not  very  small  compared  to  unity.  Therefore,  the  latter 
must  be  correct  —  most  of  the  energy  the  ions  receive  must  come  directly 
from  the  electromagnetic  fields  (induced  as  well  as  applied)  rather  than 
from  close  collisions  with  electrons.  (For  steady  fields,  the  ion  must 
receive  its  energy  from  the  electric  field  alone  because  a  steady 
magnetic  field  can  only  deflect  a  charged  particle  rather  than  impart 
energy  to  it.  Consequently,  high  specific  impulse  in  a  steady  ac¬ 
celerator  requires  high  voltage.) 

Although  the  electron  cannot  transfer  much  energy  to  an  ion 
during  an  electron-ion  collision,  the  direction  of  motion  of  the  electron 
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is  appreciably  changed  by  the  collision.  Nevertheless,  when  the  mean 

free  path  for  the  electron  is  much  greator  than  the  radius  of  gyration 

of  the  electron  in  the  magnetic  field,  (ie,  when  the  Hall  parameter 

X 

=  arfr  =  —  is  much  greater  than  one),  the  mean  motion  of  the  electrons 
6  rg 

is  nearly  the  same  as  if  they  were  collisionless.  Between  collisions 

-•  -*  p 

the  mean  velocity  of  the  electrons  is  equal  to  ExB/B  .  A  collision 

merely  shifts  the  center  of  gyration,  but  the  velocity  remains  at 

^  2 

EXB/B  .  The  shift  in  the  center  of  gyration  is  of  the  order  of  the 
gyro  radius;  so  as  long  as  the  electron  moves  much  further  than  a  gyro 
radius  between  collisions,  the  mean  motion  of  the  electrons  will  be 
nearly  the  same  as  if  the  electrons  were  coliisionless.  In  Cheng's 
experiments,  the  magnetic  field  was  of  the  order  of  one  tesla  (10^  gauss). 
The  electron  temperature  is  not  known,  but  even  if  it  were  as  low  as 
5  ev,  the  resulting  Hall  parameter  would  be  greater  than  100.  There¬ 
fore,  the  motion  of  the  electrons,  as  well  as  the  ions,  is  little  in¬ 
fluenced  by  collisions. 

Although  collisions  are  not  a  primary  mechanism  for  accelerating 

the  ions,  collisions  are  important  for  producing  the  plasma.  The 

following  more  detailed  examination  of  the  processes  that  occur  within 

the  high  specific  impulse  accelerator  indicates  how  there  may  be 

sufficient  collisions  to  produce  a  plasma  but  not  for  effective  transfer 

of  energy  between  electrons  and  ions. 

The  ion  density  throughout  most  of  the  accelerator  is  on  the  order 
15  7 

of  10  ions/ cm  ,  but  the  neutral  particle  density  near  the  point  of 

gas  injection  must  be  much  greater.  A  speculation  of  the  heavy  particle 
density  within  Cheng's  accelerator  is  given  in  figure  6.  This  distri- 
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bution  was  deduced  from  the  measured  mass  flow  rate,  the  one-dimensional 
continuity  equation,  a  measured  ion  exit  velocity,  and  a  speculated 
variation  of  axial  velocity  with  axial  position.  The  measured  velocity 
at  the  exit  is  approximately  10^  m/s  whereas  the  velocity  near  the 

region  of  gas  injection  should  be  of  the  order  of  the  thermal  velocity 

4 

of  the  gas:  less  than  10  m/s.  The  velocity  curve  shown  in  figure  6 
is  merely  a  speculation  of  how  the  velocity  increases  from  the  order  of 

4  6 

10  m/s  to  10  m/s;  the  large  difference  in  magnitude  from  the  inlet  to 

the  exit  is  the  important  feature  rather  than  the  shape  of  the  curve. 

With  a  hydrogen  flow  rate  of  the  order  of  1  g/s  and  with  a  cross  sectional 

2 

area  of  approximately  10  cm  ,  the  particle  densities  obtained  from  this 
velocity  curve  are  as  shown  on  figure  6.  The  density  near  the  inlet 
is  two  orders  of  magnitude  higher  than  the  density  throughout  the  rest 
of  the  accelerator. 

Most  of  the  ionization  should  occur  in  the  high  density  region. 
After  a  particle  is  ionized,  it  suddenly  experiences  the  high  acceler¬ 
ating  force  of  the  electromagnetic  fields.  (Note  that  the  ionization 
energy  and  the  thermal  energy  of  the  gas  at  the  inlet  are  small  compared 
to  the  kinetic  energy  of  the  ions  at  the  accelerator  exhaust).  While 
the  ion  has  low  energy,  it  should  continue  to  collide  frequently,  but 
as  it  receives  more  and  more  energy  its  collision  cross  section  decrease; 
rapidly  and  it  moves  into  a  low  density  region.  By  the  time  a  hydrogen 
ion  ha-  received  one-tenth  of  its  energy  and  has  moved  into  the  low 
density  region,  its  mean  free  path  for  collisions  with  other  ions 
(shown  in  figure  7)  is  already  several  orders  of  magnitude  larger  than 
the  size  of  the  accelerator.  Furthermore,  in  this  low  density  region, 
there  are  insufficient  electron-ion  collisions  to  transfer  an  appreciable 
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energy  between  electrons  and  ions.  Therefore,  the  hydrogen  ions 
entoring  the  low  density  region  should  become  effectively  runaway 
particles,  falling  for  all  practical  purposes  freely  through  the  electro¬ 
magnetic  fields.  (Since  they  receive  most  of  their  energy  in  this  region 
where  they  are  essentially  freely  falling,  it  is  important  that  theories 
for  these  accelerators  model  this  "free  fall  region"  correctly.) 
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Ill .  THEORETICAL  MODELING  OF  THE  ELECTRIC  PLASMA  ACCELERATORS 


The  theoretical  nodels  used  to  describe  the  electric  plasma 
accelerators  must  reflect  the  changes  in  the  accelerating  mechanisms 
described  in  tbe  previous  section. 

A  typical  conventional  arc jet  is  shown  in  figure  2.  Most  of  the 
electrical  energy  is  converted  to  thermal  energy  in  the  constricted 
part  of  the  nozzle  where  the  pressure  and  the  electrical  resistance  are 
relatively  high.  For  this  constricted  region,  one  can  model  the  energy 
conversion  and  the  gasdynamics  with  continuum  fluid  equations  with  the 
ohmic  heating  term,  J-E,  included  in  the  energy  equation  but  with  the 

—•  -4  ^ 

electromagnetic  force  terms,  pE  and  JxB,  neglected  in  the  momentum  equa¬ 
tion.  Theoretical  analyses  for  the  constricted  region  using  this  model 

5  6 

have  been  presented  by  Stine  and  Watson^  and  by  Vat son  and  Pegot  . 

There  is  relatively  little  electrical  energy  converted  to  thermal  energy 

in  the  expanded  part  of  the  nozzle,  and  this  region  may  be  approximated 

using  the  continuum  fluid  equations  with  none  of  the  electromagnetic  terms. 

A  typical  MPD  arcjet  is  shown  in  figure  8.  In  the  MPD  arc jet 

regime,  the  charged  particles  are  given  an  appreciable  deflection  between 

collisions  by  the  magnetic  field.  The  resulting  force  on  the  electron, 

expressed  by  the  J  xB  term  is  in  the  same  direction  as  the  resulting 
© 

force  on  the  ion,  expressed  by  theJ^xB  term,  and  the  otal  force, 

expressed  by  JxB,  is  not  negligible  compared  to  the  pressure  force 

as  it  was  in  the  conventional  arcjet.  The  forces  on  the  charged  par- 

tides  due  to  the  electric  field,  p  E  and  p  E,  are  opposing,  and  because 

0  1 

there  are  sufficient  collisions  such  that  the  electrons  and  ions  can 
transfer  energy,  the  net  force,  pE,  is  negligible  compared  to  the 


JxB  force  in  this  regime.  Therefore  the  theoretical  models 
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for  the  MPD  arcjet  regime  include  the  JxB  term  in  the  continuum  momentum 
equation  but  neglect  the  pE  term. 

A  number  of  attempts  have  been  made  to  obtain  solutions  for  the 

properties  of  the  MPD  arcjet  using  this  theoretical  model.  These  attempts 

7 

have  been  reviewed  recently  by  Nerheim  and  Kelly  who  concluded  that 
"the  thrust  production  mechanisms  are  poorly  understood  as  indicated  by 
the  number  of  conflicting  assumptions,  approximations,  and  processes  that 
have  been  proposed."  Although  the  only  forces  considered  for  these  attempts 
have  been  the  pressure  forces  and  the  JxB  forces,  an  accurate  solution  of 
the  pressure,  current,  and  magnetic  fields  throughout  the  accelerator  has 
not  been  obtained,  and  the  various  assumptions  and  approximations  referred 
to  by  Nerheim  and  Kelly  are  those  that  are  made  concerning  these  fields. 

In  some  analyses  the  JxB  forces  are  considered  as  constraint  forces  and 
the  propellant  is  assumed  to  be  constrained  to  flow  parallel  to  the  mag¬ 
netic  field  so  that  the  magnetic  field  surfaces  act  as  a  nozzle.  In 
other  analyses  the  JxB  forces  are  divided  into  components  and  each  compo¬ 
nent  is  treated  separately.  The  axial  component  is  called  the  "blowing 
force"  and  the  radial  component  is  called  a  "pumping  force."  When  there  is 
an  applied  magnetic  field  or  if  the  properties  are  not  axisymmetric,  there 
is  also  an  azimuthal  component.  There  is  considerable  controversy  re- 
garding  the  relative  importance  of  these  various  components  of  the  JxB 
force  and  regarding  how  they  should  be  calculated.  Nevertheless,  all  of 
the  analyses  for  the  MPD  arcjet  regime  use  this  same  model  wherein  the 

^  >4  ^ 

JxB  forces  are  included  and  the  pE  forces  are  neglected  in  the  continuum 
momentum  equation . 

In  the  direct  electromagnetic  acceleration  regime,  the  electrons 
cannot  transfer  an  appreciable  fraction  of  their  energy  to  the  ious,  and 
the  ions  receive  a  major  part  of  their  energy  from  falling  through 
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the  electric  field.  Therefore,  the  p^E  tern  becomes  the  dominant 
force  term  for  imparting  energy  to  the  ions.  Furthermore,  throughout 
most  of  the  accelerator,  neither  the  mean  free  path  of  the  ions,  shown 
in  figure  7  *8  a  function  of  ion  energy,  nor  the  ion  Larmor  radius  is 
small  compared  to  the  size  of  the  typical  accelerator;  so  the  continuum 
equations  are  inappropriate.  Therefore,  the  basic  model  used  for  the 
MPD  arcjet  regime  which  employs  the  continuum  equations  and  neglects 
acceleration  directly  due  to  the  electric  field  is  not  appropriate  in 
the  interesting  regime  of  the  direct  electromagnetic  accelerators. 

The  theoretical  model  used  herein  for  the  high  specific  impulse 
accelerators  operating  in  the  regime  of  direct  electromagnetic  accel¬ 
eration  gives  prime  consideration  to  the  "free  fall  region"  described 
in  the  previous  hhapter  because  the  ions  receive  most  <*1  their  energy 
in  this  region.  The  high  density  region  near  the  point  of  gas  in¬ 
jection  is  considered  only  as  a  source  of  plasma  for  the  upstream 
edge  of  the  "free  fall  region".  The  plasma  in  the  "free  fall  region" 
is  considered  as  collisionless  and  this  collisionless  plasma  is  sim¬ 
ulated  with  a  computer  as  described  in  the  next  chapter. 

In  1964,  in  discussions  with  Professor  R.  Jahn  of  Princeton  and 
Dr.  R.  John  of  AVCO,  the  author  proposed  that  a  regime  of  direct  electro¬ 
magnetic  acceleration  would  exist  and  suggested  a  simple  modtil  to 
describe  the  qualitative  features  of  axisymmetric  accelerators  in  this 
regime.  The  model  treats  individual  particle  trajectories  rather  than 
the  continuum  equations,  and  the  ions  are  assumed  to  receive  their 
energy  directly  from  the  electric  fields.  At  about  the  same  time, 

Q 

Stratton0  formulated  a  very  similar  model  that  included  two  additional 
assumptions  regarding  the  form  of  the  magnetic  field  and  the  size  of  the 
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ion  Larmor  radius.  These  additional  assumptions  allowed  him  to  make 
quantitative  predictions  for  the  accelerators.  These  two  models  will 
be  described  and  compared  with  the  results  of  the  computer  simulation 


after  a  discussion  of  the  computer  simulation  study  to  be  described 


IV.  COMPUTER  SIMULATION  OF  THE  DIRECT  ELECTROMAGNETIC  ACCELERATOR 


A.  Introduction 

In  order  to  obtain  a  solution  for  the  notion  of  the  plasma  in  the 
"free  fall  region"  of  the  high  specific  Impulse  accelerator^  a  computer 
was  used  to  simulate  the  plasma.  Tills  simulation  was  achieved  by  rep¬ 
resenting  the  many  ions  and  electrons  within  the  accelerator  by 
several  thousand  particles  and  by  following  the  motion  of  these  particles 
in  time;  i.e.,  by  the  Lagrangian  technique.  Time  was  divided  into  small 
discrete  increments  and  the  following  computation  procedure  was  used  to 
provide  a  time  history  of  the  fields  and  particle  motions. 
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REPEAT  FOR  EACH  INCRBCENT  IN  TIME 


A  lucid  description  of  the  Lagrangian  technique  vas  given  by 
Buneman  and  Dunn  in  reference  9«  The  technique  has  been  previously 
employed  primarily  to  study  electrostatic  waves  in  a  one-  or  two-di¬ 
mensional  ,  low  density,  boundless  plasma  or  to  study  the  ion  trajectories 
and  electron  neutralization  within  ion  guns .  In  both  applications,  the 
self- induced  magnetic  field  has  been  neglected.  A  survey  of  the  ap¬ 
plications  of  this  technique  is  given  in  reference  10.  The  computer 
simulation  presented  herein  differs  from  the  majority  of  previous 
applications  in  that  boundaries  and  self-induced  magnetic  fields  are 
included,  the  domain  is  an  axlsymmetric  space  rather  than  a  one-  or 
two-dimensional  cartesian  space,  and  the  simulated  plasma  densities  are 
much  greater. 

B.  Accelerator  Configurations,  Boundary,  and  Initial  Conditions 

The  two  accelerator  configurations  used  in  the  computer  simu¬ 
lation  are  shown  in  figure  9  and  were  selected  to  simulate  the  acceler- 

2  U 

ators  being  tested  by  Cheng  (shown  in  figure  3)  and  by  Clark  and  Jahn 

(shown  in  figure  4).  Since  detailed  experimental  measurements  by  Clark 
and  Jahn  showed  that  the  plasma  properties  are  axi symmetric,  axisymmetry 
was  assumed  for  the  computer  simulation. 

The  insulator  extending  downstream  of  the  anode  on  the  models 
for  the  computer  simulation  were  added  to  simplify  tho  mesh  configura¬ 
tions  (shown  in  figure  10)  for  the  computer  solution  of  the  electro¬ 
magnetic  fields.  The  insulator  could  be  moved  out  to  a  larger  diameter 
with  new  mesh  points  added  to  fill  the  added  space.  However,  the 
program  presented  uses  all  of  the  word  memory  of  the  computer  which  was 
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employed;  so  if  more  mesh  points  are  used,  the  number  of  ions  and 
electrons  that  simulate  the  plasma  would  have  to  be  reduced  or  &  larger 
computer  employed.  Sven  though  the  added  insulator  causes  some  reduction 
in  perforaance,  the  solutions  are  adequate  to  display  the  significant 
qualitative  features  of  the  accelerator  and  to  check  the  computer 
simulation  techniques. 

It  was  assumed  that  a  constant  potential  was  applied  between  the 
electrodes  and  that  the  power  supply  was  electrically  floating.  In  the 
computer  simulation  these  two  conditions  were  used  to  calculate  the 
surface  charges  on  the  electrodes  at  each  time  step.  The  condition 
that  the  power  supply  be  electrically  floating  combined  with  the  con¬ 
servation  of  charges  gives  the  requirement  that  the  net  transfer  of 
electrons  from  the  external  circuit  to  the  plasma  (via  the  electrodes) 
must  be  the  same  as  the  net  change  in  the  total  charge  of  the  external 
circuit.  In  the  computer  simulation,  the  sum  of  the  electrode  surface 
charges  was  required  to  equal  the  net  transfer  of  the  electrons 
from  the  electrodes  to  the  plasma.  The  condition  that  the  potential 
difference  between  the  electrodes  be  constant  gives  the  requirement 
for  tho  difference  between  the  surface  charges  of  the  electrodes.  These 
two  requirements  are  sufficient  to  uniquely  determine  the  electrode 
surface  charges.  Note  that  the  floating  power  supply  condition  auto¬ 
matically  assures  that  for  steady  state  operation,  the  plasma  leaving 
the  accelerator  will  be  neutral.  If  the  electron  flow  rate  from  the 
circuit-accelerator  combination  is  greater  then  the  ion  flow  rate,  then 
an  excess  of  ions  builds  up  in  the  circuit-accelerator  combination 
increasing  the  electric  potential  until  the  electron  and  ion  flow  rates 
become  the  same. 
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No  magnetic  field  was  applied  —  the  only  magnetic  field  It  1 

that  induced  by  the  particle  notions.  Since  the  model  is  axisymaetrlc  | 

with  no  azimuthal  current,  the  induced  magnetic  field  is  in  the 

I 

azimuthal  direction  only. 

I 

A  neutral  gas  mas  assumed  to  be  injected  through  the  insulator 

vails  surrounding  the  cathode  and  the  gas  was  arbitrarily  assumed  to 

become  ionized  in  the  regions  shown  in  figure  9.  (A  simulation  for 

Clark  and  Jahn's  accelerator  where  the  ionization  region  was  assumed  to 

be  at  the  anode  surface  -  considerably  different  than  the  region  shown 

in  figure  9  ~  illustrated  that  the  location  of  the  ionization  region 

influences  the  efficiency  Of  the  accelerator  somewhat,  as  described  in 

a  later  section  entitled  "Discussion",  but  it  did  not  change  the 

qualitative  features  described  in  a  later  section  entitled  "Primary 

Results  of  the  Computer  Simulation").  The  neutral  particles  were  not 

included  in  the  simulation  —  rather  electron-ion  pairs  were  created 

in  the  assumed  'onization  region  at  a  rate  corresponding  to  the  mass 

flow  rate  assuming  100$  ionization  of  the  gas.  (Note  that  the  gas 

can  be  assumed  to  be  completely  ionized  even  though  the  hydrogen  ions 

are  essentially  collisionless  in  the  low  density  region  because  the 

cross  section  for  ionization  of  the  hydrogen  atom  is  much  larger  than 

the  energy  exchange  cross  sections  for  high  energy  ion-ion  and  electron- 

ion  encounters.  For  example,  the  cross  section  for  ionization  of  the 

11  -17  2 

hydrogen  atom  given  by  Jahn  is  greater  than  10  cm  whereas  the 

4 

ion-ion  cross  sections  for  10  ev  ions  —  the  energy  corresponding  to 

c  ^2^)  O 

a  specific  impulse  of  10?  seconds  —  is  approximately  10  cm  . )  The 
simulations  were  found  to  be  insensitive  to  the  assumed  thermal  velocity 
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of  the  newly  created  ions  and  electrons  as  long  as  the  assumed  thermal 
energy  was  much  less  than  the  energy  of  the  ions  at  the  accelerator 
exhaust . 

The  current  drawn  from  the  electrode  was  assumed  to  be  proportional 
to  the  axial  electric  field  at  the  cathode,  and  a  high  axial  field  was 
assumed  to  induce  a  very  large  current  to  simulate  a  high  electrical 
conductivity  near  the  cathode.  The  currents  did  not  become  excessive 
because  the  axial  electric  fields  became  small  as  described  in  a  later 
section  entitled  "Discussion". 

When  the  particles  hit  the  boundaries,  the  encounters  can  be 
treated  arbitrarily  as  specular,  diffuse,  or  any  combination  of  these. 

For  the  calculations  to  be  presented,  the  particles  were  arbitrarily  as¬ 
sumed  to  lose  half  of  their  energy  with  the  angle  of  reflection  equal 
to  the  angle  of  incidence.  The  exceptions  are  that  the  electrons  were 
absorbed  at  the  anode  and  the  ions  were  considered  to  recombine  when 
striking  the  cathode.  Except  for  the  wall  encounters,  the  plasma  was 
assumed  to  be  collisionless;  no  hard  collisions  between  particles  were 
included . 

The  simulation  was  started  with  the  accelerator  void  of  particles 
and  with  the  simultaneous  initiation  of  mass  flow  and  current.  The 
computer  simulation  passed  through  a  transient  phase  and  then  came  to 
a  steady  state  operation.  The  simulation  was  formulated  to  represent 
the  steady  state  phase;  the  transient  phase  was  not  simulated  accurately — 
rather  the  simulation  was  started  with  an  initial  condition  and  allowed 
to  progress  in  time  merely  as  a  convenient  means  for  obtaining  the 
steady  state  operation.  The  modification  of  the  program  that  would 
be  required  to  simulate  the  transient  phase  properly  is  described  in 
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a  later  section  entitled  "Details  of  the  Simulation  Technique", 

C.  Equations  of  Motion 

The  motion  of  each  particle  is  calculated  using  approximations  to 
the  following  exact  equations  of  motion. 

e,-»  -•  «*. 

dv  =  Z-(E  +  v  X  B)  dt 
m'  ' 

^  — • 

dx  =  v  dt 

where  v  =  velocity,  m/s 

Z  =  number  of  charges  on  the  particle 
(  =  -1  for  electrons) 

e  =  charge  of  singly  ionized  ion,  coulombs 
m  =  mass  of  the  particle,  kg 
E  =  electric  field,  V/m 
B  =  magnetic  field,  tesla 
t  =  time,  s 

x  =  position  vector,  m 
The  approximations  employed  are: 

(1)  The  differential  equations  are  approximated  by  finite  dif¬ 

ference  equations. 

(2)  The  azimuthal  velocity  is  neglected.  (The  velocities  induced 

by  electromagnetic  acceleration  art  much  greater  than  the 
thermal  velocities  of  the  injected  particles,  and  these 
induced  velocities  are  in  the  axial  tmd  radial  direction 
only) 
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With  these  approximations,  the  equations  of  motion  as  derived  by 

_  ,  12  , 

Hockney  become: 
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where  v  =  axial  velocity  at  the  new  time,  m/s 
2 

new 

v  =  axial  velocity  at  the  previous  time,  m/s 
Zold 

v  ss  radial  velocity  at  the  new  time,  m/s 

new 

v  =  radial  velocity  at  the  previous  time,  m/s 
rold 

r  =  radial  position  at  the  new  time,  m 
new  ’ 

=  radial  position  at  the  previous  time,  m 
z  =  axial  position  at  the  new  time,  m 

zq1^  =  axial  position  at  the  previous  time,  m 

At  =  the  time  increment, s 

T)  =  At  Z— 

m 

e  =  t$/2 


These  finite  difference  equations  were  formulated  to  yield  exact 
solutions  for  two-dimensional  motion  when  the  electric  end  magnetic 
field  are  constant,  and  they  are  good  approximations  as  long  as  the  fields 
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have  relatively  saall  variations  during  the  tiae  increment. 


D.  Field  Equations 

The  electric  and  nagnetic  fields  used  in  the  equations  of  action 
given  in  the  previous  section  were  obtained  using  approximations  to  the 
following  exact  field  equations. 


E  =  -•$  -  dA/dt 


B  =  V  X  A 

* =  fdk dVolume 

all  space 

*  =  /ini dVolUBe 

all  space 
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electric  field,  V/m 
electric  potential,  V 
magnetic  potential  rector.  Weber/m 
magnetic  field,  tesla 

3 

charge  density,  coulombs/m 

2 

current  density.  A/m 

electric  permittivity,  farad/m 

2 

magnetic  permeability,  N/A 
distance  between  the  point  of  measure¬ 
ment  of  $  or  A  and  the  location  of 
the  charge  or  current  element 


The  approximations  employed  for  the  field  equations  are: 


(l)  The  plasma  and  field  properties  are  assumed  to  be  axlsymmefrlc 


r  4 

(The  experimental  measurements  by  Clark  and  Jahn  shoved 
that  the  properties  mere  axlsymmetrlc  in  their  accelerator) 

(2)  The  plasma  and  field  properties  are  assumed  to  attain  a 

steady  state.  (The  experimental  measurements  by  Clark  and 

4 

Jahn  shoved  that  the  properties  attained  a  steady  state) 

(3)  The  fields  are  assumed  to  be  smooth  enough  to  be  approximated 

by  linear  interpolation  of  field  values  calculated  at 
relatively  few  mesh  points  within  the  accelerator. 

Using  the  first  approximation,  one  can  obtain  the  magnetic  field  from 
Ampere's  circuital  lav: 

Bg=  pl/2nr 

where  r  -  radial  distance  from  the  axis  of  symmetry,  m 

I  =  the  current  that  passes  through  the  circle  formed 
by  rotating  the  point  of  measurement  about 
the  axis  of  symmetry,  A 

From  the  second  approximation,  the  electric  field  is  given  by: 


E  =  -V$ 


where  the  electric  potential,  using  the  third  approximation, 
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by: 


$  = 
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mn  n 


n=l 


is  given 


where  $  =  the  electric  potential  at  the  m  mesh  point 

m 

th 

C  =  the  electric  potential  at  the  m  mesh  point  induced 
mn 

th 

by  a  unit  charge  in  the  n  mesh  cell 

th 

pQ  s  the  charge  density  in  the  n  mesh  cell 
N  =  the  total  number  of  mesh  cells. 


The  aesh  configurations  are  shown  in  figure  10.  Note  that  the  aesh 
points  are  actually  circles  around  the  axis  of  symmetry  and  the  aesh 
cells  are  actually  toroids  around  the  axis  of  syaaetry.  The  detailed 
technique  for  calculating  the  electric  potentials  at  the  aesh  points 
is  described  in  the  next  section. 

E.  Details  of  the  Simulation  Technique 

The  computer  simulations  were  obtained  using  a  modest  size 
computer  —  an  IBBt  709^  with  a  32K  \iot d  direct  access  core  —  with  the 
computer  program  given  in  the  appendix.  The  primary  purpose  of  this 
work  is  to  introduce  and  Illustrate  simulation  concepts  and  to  check 
the  qualitative  features  of  the  simple  models  proposed  by  Stratton  and 
the  author  rather  than  to  present  highly  detailed  solutions.  Neverthe¬ 
less,  the  computer  program  presented  in  the  appendix  can  easily  be 
modified  to  give  additional  detail  and  more  accuracy  with  the  larger 
computers  and  with  more  computation  time. 

The  simulation  techniques  employed  in  the  computer  program  that 
have  been  previously  employed  by  others  and  that  are  explained  by 

Q 

Buneman  and  Dunn  are  as  follows.  The  many  electrons  and  ions  of  the 
plasma  are  represented  by  a  far  fewer  number  of  representative  particles 
—  up  to  2000  electrons  and  2000  ions  for  the  programs  presented  in  the 
appendix.  The  electric  and  magnetic  fields  were  calculated  from  the 
particle  distributions  only  at  discrete  points  within  the  domain  —  on 
the  mesh  points  shown  in  figure  10  for  the  programs  presented.  The 
fields  between  these  points  were  obtained  by  linear  interpolation.  A 
highly  reduced  ion-electron  mass  ratio  was  used  —  for  the  calculations 
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presented  herein,  usually  a  ratio  of  10  to  1,  These  techniques  were 
checked  for  this  problem  as  follows.  Tests  wherein  the  number  of  par¬ 
ticles  used  to  represent  the  plasma  was  varied  indicated  that  several 
thousand  particles  is  an  adequate  statistical  sample  to  use  for  this 
two-dimensional  problem.  Tests  with  the  mesh  size  reduced  to  half  the 
size  shown  in  figure  10  indicated  that  the  course  mesh  was  adequate 
because  the  fields  were  smoothly  varying  for  this  problem.  Tests  with 
mass  ratios  varying  from  2  to  32  Indicated  that  the  steady  state 
solutions  were  only  slightly  sensitive  to  a  variation  of  the  mass  ratio 
from  2  to  8  and  Insensitive  to  the  variation  of  mass  ratio  from  8  to 
32. 

The  reduced  mass  ratio  is  used  in  the  computer  simulation  to 
reduce  the  computation  time.*  The  time  step  used  in  the  simulation 
must  be  smaller  than  the  inverse  of  the  plasma  frequency  in  order  that 
the  numerical  solution  be  stable. 


At  <  - 

CD 


where  At  =  the  time  step  used  in  the  simulation,  s 

n  e2  \  1/2 


0)p  =  the  plasma  frequency 


)  ,  radians/ s 

e  m  /  ’ 

o  e/ 

nfi  =  electron  number  density,  electrons/m-3 


-e 


e  =  charge  of  a  singly  ionized  ion,  coulombs 

C  =  electric  permittivity  of  free  space,  farads/m 
o 


m  =  electron  mass,  kg 
© 


(Note  that  the  plasma  frequency  is  a  natural  oscillation  frequency  of 
the  electrons,  analogous  to  the  natural  oscillation  frequency  of  a  spring' 
mass  system;  hence,  one  can  see  the  physical  reason  why  the  computations 


become  unstable  if  the  time  step  exceeds  the  inverse  of  the  plasma 
frequency) .  By  increasing  the  electron  mass  to  one-tenth  the  ion  mass, 
the  time  step  can  be  increased  and  the  computation  time  reduced  more 
than  one  order  of  magnitude.  (Note  that  the  increased  electron  mass 
changes  the  plasma  frequency;  so  if  one  wanted  to  simulate  the  electron 
fluctuatiois  in  an  unsteady  plasma,  one  would  have  to  use  the  proper 
electron  mass.) 

A  new  technique  for  simulating  the  steady  flow  of  high  density 
plasmas  was  examined  and  employed  for  the  simulations  presented  herein. 
Since  the  minimum  time  step  is  inversely  proportional  to  the  square  root 
of  the  electron  number  density,  zs  shown  above,  the  computation  time 
is  proportional  to  the  square  root  of  the  electron  number  density; 
therefore,  the  simulation  of  dense  plasmas  using  the  Lagrangian  technique 
described  above  had  not  been  attempted  previously.  However,  the  sim¬ 
ulations  presented  herein  for  the  steady  dense  plasma  flows  within 
the  accelerators  were  obtained  with  reasonable  computation  times 
by  artificially  increasing  the  electric  permittivity,  which  reduces 
the  plasma  frequency  and  reduces  the  computation  time.  (Note  that  this 
technique  Is  similar  to  artificially  increasing  the  electron  mass  in 
that  both  artificially  reduce  the  plasma  frequency.  Neither  technique 
can  be  used  for  simulating  plasmas  where  the  electron  fluctuations  must 
be  considered) . 

Some  understanding  of  the  effect  of  artificially  increasing  the 
electric  permittivity  can  be  obtained  by  considering  the  flow  of  an 
inviscid  plasma  through  an  insulated  channel.  When  an  electric  field 
is  applied  transverse  to  the  flow,  the  electrons  and  ions  are  deflected 
in  opposite  directions  toward  the  channel  walls  until  plasma  sheaths 
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are  formed  which  counterbalance  the  applied  electric  field.  The 
particle  trajectories  are  not  altered  greatly  by  the  applied  electric 
field  if  the  relative  shift  of  electrons  and  ions  to  form  the  plasma 
sheaths  is  small.  The  distance  the  electrons  and  ions  shift  to  counter¬ 
balance  the  applied  field  is 


where  6  =  the  distance  the  electrons  and  ions  shift 

toward  the  walls  to  form  the  plasma  sheath,  m 
=  electric  permittivity  of  free  space,  farads/m 
e  =  charge  of  a  singly  ionized  ion,  coulombs 
♦  s  the  electric  potential  applied  across  the 
channel,  V 

n  =  electron  number  density,  electrons/mJ 
e 


The  deflection  of  the  individual  particles  is  very  small  compared  to 
the  channel  width,  W,  when 


Note  that  this  condition  is  the  same  as  the  condition  that  the  applied 
potential  be  much  smaller  than  the  potential  that  would  be  produced 
by  removing  all  of  the  ions  from  the  channel,  or,  equivalently,  that 

A  IPi-Pe1 

the  relative  charge  density,  p  - - ,  be  much  less  than  unity. 

Pr 

Artificially  increasing  eQ  pennies  the  electrons  and  ions  to  shift  more, 
creating  an  artificially  large  relative  charge  density.  However,  little 
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error  in  the  particle  trajectories  is  produced  by  increasing  e  as 

long  as  the  above  condition  ia  net;  i.e.,  as  long  as  the  relative 

charge  densities  are  small  compared  to  unity. 

The  ion  and  electron  number  densities  within  the  accelerators 

15  a 

simulated  are  of  the  order  of  10  particles/cm.  However,  the  net 

charge  density,  the  difference  between  the  ion  and  electron  number 

q  v 

densities,  is  of  the  order  of  10  charges/ cm  .  Therefore;  the  relative 

charge  density,  3,  is  of  the  order  of  10~^.  Although  a  frequently 

used  technique  for  treating  plasmas  with  smell  relative  charge  densities 

is  to  artificially  set  $  identically  equal  to  zero,  assuming  perfect 

neutrality,  the  technique  employed  herein  was  to  artificially  increase 

the  electric  permittivity,  e permitting  jj  to  increase  from  the  order 

of  10  ^  to  the  order  of  10  1.  Tests  that  were  performed  varying  the 

electric  permittivity  illustrated  that  the  ion  and  electron  distributions 

adjusted  slightly  to  give  very  nearly  the  same  electric  potential  despite 

the  variation  in  electric  permittivity.  (Only  slight  adjustments  of 

the  particle  trajectories  were  required  to  give  these  adjustments  in 

electron  and  ion  distributions.)  In  fact,  the  qualitative  features 

of  the  simulation  were  satisfactorily  given  even  with  the  relative 

charge  densities  as  high  as  1/3. 

13 

Butler,  Cook,  and  Morse  have  assumed  complete  neutrality, 

$  =  0,  in  order  to  obtain  a  computer  solution  of  a  coaxial  accelerator 
such  as  Cheng's.  Whenever  complete  neutrality  is  assumed,  the 
electric  field  cannot  be  calculated  from  the  charge  distributions  — 
rather,  it  must  be  calculated  from  the  motion  of  the  plasma.  Butler, 
Cook,  and  Morse  calculated  the  particle  motion  by  assuming  the  plasma 
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was  frozen  to  the  magnetic  field  lines  and  by  treating  the  magnetic 
pressure  as  the  dominant  accelerating  force.  (The  electric  fields  are 
not  even  used  in  this  type  of  calculation  even  though  in  actuality  the 
ions  receive  all  of  their  energy  from  the  electric  fields.)  Their 
technique  should  take  less  computing  time  for  simulating  plasmas  that 
are  frozen  to  field  lines.  However,  it  is  shown  in  a  later  section 
entitled  "Primary  Results  of  the  Computer  Simulation"  that  the  plasma 
is  not  frozen  to  the  magnetic  field  lines  within  Cheng's  accelerator. 

The  potential  field  at  each  mesh  point  was  calculated  by  super¬ 
imposing  the  potentials  induced  at  the  mesh  point  by  the  surface 

charges  and  by  the  charges  within  all  of  the  mesh  cells  rather  than 

12  14 

calculated  from  the  technique  developed  by  Hockney  ’  .  (The  charges 

were  assigned  to  the  mesh  cells  using  the  method  described  by  Birdsall 
and  Fuss  J .)  The  potential  induced  on  all  mesh  points  by  a  unit  charge 
in  one  mesh  cell  or  by  a  unit  surface  charge  at  an  electrode  was 
calculated  once  at  the  start  of  the  program,  and  these  inverse  capaci¬ 
tances  were  merely  multiplied  by  the  actual  charge  in  each  mesh  cell  or 
the  actual  surface  charge  to  calculate  the  electric  potentials  at  each 
time  step.  This  technique  is  particularly  suitable  for  problems  with 
mixed  boundary  conditions,  such  as  for  intermittent  conductors  and  in¬ 
sulators,  and  for  problems  with  less  than  1000  mesh  points.  For 
problems  with  more  mesh  points  and  with  the  special  boundary  conditions 
cited  by  Hockney,  the  technique  he  presented  would  take  less  computer 
time.  The  solution  for  the  potential  fields  using  the  technique  just 
described  for  the  mesh  configuration  shown  in  figure  10  required  5 
milliseconds  on  an  IBM  709^* 
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The  steady  magnetic  field  was  easily  determined  for  this 


The  current  density  was  calculated  at  each  mesh  point  from  the 
particle  velocities  for  the  particles  inside  the  mesh  cell. 

To  calculate  unsteady  fields^  one  could  use  the  magnetic  poten¬ 
tial  vector  with  the  equation 


and  the  same  program  for  the  solution  of  Poissons  Equation  that  is  em¬ 
ployed  to  solve  for  the  electric  potentials  could  be  used  for  the 
solution  of  the  magnetic  potentials.  Such  calculations  as  these  would 
be  necessary  in  order  to  properly  simulate  the  transient  phase  of  the 
accelerator  operation. 

The  computation  times  required  to  achieve  steady  state  solutions 
for  the  plasma  accelerators  were  approximately  oO  minutes  on  the  IBM 
7094  computer  using  several  thousand  representative  particles  and  using 
the  mesh  configuration  shown  in  figure  10. 

F.  Display  Techniques 

A  computer  program  was  developed  to  display  the  calculated  particle 
motions,  the  electric  pc  .ential,  and  the  magnetic  field  simultaneously 
on  a  cathode  ray  tube.  The  cathode  ray  tube  display  was  recorded  on 
motion  pictures;  a  typical  frame  of  the  motion  pictures  is  shown  in 

figure  11.  Since  the  model  is  assumed  axisymmetric,  the  properties  on 
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any  plane  through  the  axis  are  the  same.  The  motion  of  the  particles 
is  shown  on  one  such  plane  in  the  center  figure,  which  is  schematically 
shown  as  a  cross  section  of  the  accelerator.  The  same  plane  has  been 
translated  to  the  upper  left  to  show  the  electric  potential  and  to 
the  lower  right  to  show  the  magnetic  field.  Since  no  magnetic  field 
is  applied,  the  model  is  axisymmetric,  and  the  azimuthal  velocities  are 
assumed  to  be  zero,  the  only  component  of  the  magnetic  field  is 
azimuthal;  the  absolute  magnitude  of  this  azimuthal  field  is  the 
value  plotted.  The  motion  pictures  permit  one  to  observe  the  interaction 
between  the  particles  and  the  fields.  One  can  see  how  the  fields  ac¬ 
celerate  the  particles,  hov  the  particle  motion  changes  the  charge 
density  and  current  distributions,  and  how  the  changing  distributions 
modify  the  fields  that  accelerate  the  particles. 

The  calculations  for  the  simulation  and  the  display  were  not  done 
simultaneously.  The  particle  positions  and  the  electromagnetic  fields 
were  calculated  on  the  IBM  709^  and  these  quantities  were  stored  on 
magnetic  tape.  The  information  on  the  magnetic  tape  was  then  processed 
and  plotted  with  a  CDC  3200  computer-CDC  28 0  cathode  ray  tube  com¬ 
bination.  T1  is  computer-cathode  ray  tube  combination  produced  the  frames 
illustrated  on  figure  11  at  approximately  2  frames/sec.  The  computation 
time  required  to  plot  the  simulation  was  approximately  1/10  the  time  to 
generate  the  information  with  the  IBM  JOjb. 

Because  of  the  large  volume  of  information  processed,  special 
care  was  taken  to  minimize  tape  read  and  write  times,  and  parallel 
operations  were  performed  where  possible.  For  example,  the  floating 
numbers  generated  for  the  particle  positions  and  for  the  electro¬ 
magnetic  fields  were  converted  to  12  bit  integers,  and  three  such 
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integers  were  packed  into  each  36  bit  IBM  7094  word  prior  to  recording 
the  information  onto  magnetic  tape.  The  CDC  3200  was  programmed  so 
that  it  plotted  the  information  for  one  frame,  processed  the  information 
of  a  second  frame,  and  read  the  information  for  a  third  frame  simul¬ 
taneously. 

G.  Primary  Results  of  the  Computer  Simulation 

The  more  important  qualitative  features  observed  in  the  computer 
simulation  movies  are  as  follows.  Despite  the  difference  in  shape 
between  Cheng's  accelerator  and  Clark  and  Jahn’s  accelerator,  the  same 
features  are  noted  for  both.  The  self-induced  magnetic  field  causes 
a  large  fraction  of  the  electrons  emitted  from  the  cathode  to  be  de¬ 
flected  from  a  straight  path  to  the  anode  toward  a  path  parallel  to  the 
axis  of  the  accelerator.  The  increase  in  electron  density  near  the 
axis  in  turn  lowers  the  electric  potential  along  the  axis  such  that  the 
electric  field  in  the  region  of  ionization  becomes  mainly  radial  as 
shown  in  figure  12.  A  large  fraction  of  the  ions  accelerated  toward 
the  axis  by  the  radial  field  are  deflected  by  the  induced  magnetic  field 
toward  the  accelerator  exhaust  along  a  path  similar  to  the  one  shown 
in  figure  13. 

Most  of  the  electrons  that  reach  the  anode  are  the  electrons 
released  by  the  ionization  of  the  gas,  rather  than  electrons  that 
originate  at  the  cathode.  The  majority  of  the  electrons  that  originate 
at  the  cathode  leave  the  accelerator  with  the  ions  that  are  deflected 
to  the  accelerator  exhaust  forming  the  neutral  plasma  stream  that  leaves 
the  accelerator.  The  radial  current  is  not  everywhere  mainly  electron 
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current  as  in  the  MPD  arjet  —  rather  an  appreciable  part  of  the 
radial  current  is  carried  by  ions.  With  the  ionization  regions  as 
shown  in  figure  9>  approximately  half  of  the  total  radial  current  is 
ion  current . 

The  simulation  shows  that,  for  this  problem,  two  frequently  used 
concepts  must  be  modified.  The  concepts  are  (1)  that  the  p^E  force 

—4 

cancels  the  pgE  force  for  a  neutral  plasma,  and  (2)  that  plasmas  with 
infinite  conductivity  are  frozen  to  magnetic  fiexd  lines.  Although 

—4 

the  plasma  is  nearly  neutral,  the  p^E  force  is  not  strictly  cancelled 
by  the  p^E  force  as  in  a  highly  colllsional  plasma  because  in  the  direct 
electromagnetic  accelerator  the  electron  and  ion  gases  can  easily  slip 
through  one  another.  Rather,  the  p^E  force  is  balanced  partially  by 
inertia  forces  and  partially  by  the  v^  X  B  force.  Figure  14  shows  an 
ion  trajectory  for  an  ion  starting  from  rest  in  a  uniform  electric 
field  and  a  decaying  magnetic  field  —  this  trajectory  is  similar  to 
the  ion  trajectories  observed  in  the  numerical  simulation.  As  the  ion 

—4 

begins  its  motion,  the  p^E  force  is  balanced  almost  entirely  by  the 

inertia  force  whereas  during  the  last  part  of  the  trajectory  shown  in 

-»  —  — 

figure  14,  the  p^E  force  is  balanced  by  the  v.  X  B  force.  The  peE 

•4 

force  on  the  electrons  is  balanced  almost  entirely  by  the  ve  X  3  force 
since  the  inertia  of  an  electron  is  so  much  smaller  than  the  inertia  of 
an  ion.  The  typical  electron  trajectory  for  the  same  fields  as  for  the 
ion  is  also  shown  on  figure  l4. 

Although  the  electrical  conductivity  in  the  simulation  is  effec¬ 
tively  infinite  because  there  are  no  particle-particle  collisions  in¬ 
cluded  in  the  simulation,  the  fluid  is  not  frozen  to  the  magnetic  field. 
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With  the  reference  frane  fixed  in  the  accelerator,  the  magnetic  field 
lines  are  stationary  and  the  particles  are  accelerated  through  the*  on 
trajectories  such  as  shown  ir.  figure  13.  Furthermore,  since  the  velocity 
of  the  ions  and  electrons  vary  in  direction  and  magnitude  throughout  the 
accelerator,  there  is  no  transformation  of  reference  frames  possible 
where  the  electric  field  could  be  transformed  to  a  moving  magnetic  field 
in  which  the  fluid  would  be  frozen  to  the  magnetic  field  lines.  One  can 
deduce  from  the  simulation  that  it  is  the  inertia  terms  that  cause  the 
slippage  of  plasma  through  the  field  lines.  Particles  with  finite  inertia 
may  rotate  around  the  field  lines  with  a  finite  gyro  radius.  In  the  com¬ 
puter  simulation,  the  gyro  radius  for  the  ions  is  approximately  equal  to 
the  size  of  the  accelerator  and  to  the  size  of  the  volume  that  contains 
a  large  magnetic  field. 

H.  Comparison  of  the  Computer  Simulation  with  the  Simple  Models 

The  essential  features  that  were  incorporated  as  a  priori  assump- 

8 

tions  in  the  simple  models  proposed  by  Stratton  and  the  author  were  cor¬ 
roborated  by  the  computer  simulation.  These  assumed  features  are:  (l)  the 
self-induced  magnetic  field  channels  the  electrons  from  the  cathode  along 
the  axis,  causing  a  dip  in  the  potential  field  along  the  axis  —  the  ions 
then  receive  their  energy  from  the  resulting  predominately  radial  field; 
and  (2)  the  self-induced  magnetic  field  converts  the  radial  velocity  ob¬ 
tained  by  the  ions  from  the  electric  field  to  an  axial  velocity  as  shown  in 
figure  13.  These  features  are  not  assumed  for  the  computer  simulation  — 
rather,  the  majc  assumptions  for  the  computer  simulation  are  that  the 
dominant  acceleration  mechanism  can  be  simulated  with  a  collisionless 
model  and  that  sufficient  particles  are  treated  in  the  simulation  to 
comprise  a  reasonable  statistical  sample.  Because  the  two  features 
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are  found  as  a  result  of  calculating  the  self-consistent  fields  and 
particle  Motions  in  the  computer  simulation,  the  simulation  corroborates 
the  a  priori  assumptions  of  the  simple  models. 

These  features,  when  stated  in  equation  form  are: 

m  v.^ 

"V"  =  ciZe$  (0 

the  ions  can  receive  energy  only  from  the  electric  field  where 
m^  =  mass  of  the  ion  ,  kg 

v  =  the  exit  velocity  of  the  ion,  m/s 

=  a  proportionality  constant  less  than  1 

§  =  the  potential  applied  to  vhe  accelerator,  V 


Z  =  number  of  charges  on  each  ion 


e  =  charge  of  a  singly  ionized  ion,  coulombs 


axial 


ion  path 


2£iS!_  dt  =  f  2®>dr  (2) 

*1  /  "  1 


ion  path 


the  radial  velocity  received  from  the  electric  field  is  converted  to 
axial  velocity  by  the  magnetic  field  where 

Bg  =  aximuthal  magnetic  field,  tesla 
t  =  t  iir.e  ,  s 

r  =  radial  position  within  the  accelerator  ,m 

g 

Stratton  and  the  author  further  proposed  that,  for  efficient 
performance,  the  number  of  neutral  particles  injected  into  the  accelerator 
times  the  average  number  of  charges  per  ion  should  be  approximately  the 
same  as  the  number  of  electrons  emitted  at  the  cathode.  This  hypothesis  was 


not  studied  with  the  computer  simulation;  however  the  following 
argument  is  presented  to  support  its  validity:  for  steady  operation 
the  number  of  electrons  leaving  the  accelerator  must  equal  the  number 
of  positive  charges  (number  of  ions  times  the  average  number  of  charges 
per  ion)  leaving  the  accelerator.  Forcing  the  number  of  electrons 
emitted  at  tne  cathode  to  be  much  larger  would  merely  force  many  of  the 
electrons  across  the  potential  to  the  anode,  and  since  the  electrons 
cannot  effectively  transfer  their  energy  to  the  ions,  this  would  merely 
result  in  additional  heating  to  the  anode. 

Stated  in  equation  form  this  is  written 


where 

m  =  mass  flow  rate  of  injected  gas  , kg/s 
I  =  current  ,A 

C2  =  proportionality  constant  (approximately  equal  to  unity) 

If  the  relationship  between  the  current  and  the  magnetic  field  is  known 
or  assumed,  then  equation  (2)  becomes 


Vj_  =  C3I  (4) 

where  =  proportionality  constant. 

In  order  to  predict  quantitative  characteristics  for  direct 
electromagnetic  accelerators,  the  proportionality  constants  C^,  C,)}  and 
must  be  known  or  assumed.  Since  the  potential  the  ions  fall  through 
is  of  the  same  order  as  the  applied  potential,  must  be  of  the  order 
of  unity.  C must  be  of  the  order  of  unity  for  efficient  operation, 
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and  was  set  equal  to  one  for  the  computer  simulation  by  setting  the 
rate  of  ionization  equal  to  the  rate  of  electron  Emission  from  the 
cathode.  In  order  to  obtain  Stratton  assumed  the  magnetic  field 
to  be 


* 


hI_ 

2nr 


(5) 


where  p  =  magnetic  permeability,  and  assumed  the  ion  gyro  radius  is 
everywhere  equal  to  the  geometrical  radius  from  the  axis  of  the 
accelerator.  With  these  assumptions,  Stratton  obtained  a  value  of  20 
for  C^.  On  the  other  hand,  the  computer  simulation  yields  a  value  of 
of  approximately  15  for  the  model  of  Clark  and  Jahn' s  accelerator  and 
approximately  6  for  the  model  of  Cheng's  acceleratoi .  The  magnetic 
field  obtained  from  the  computer  simulation,  shown  in  figure  11,  does 
not  compare  well  with  Stratton's  assumption  given  by  liquation  5,  and 
the  typical  ion  path  shown  in  figure  13  does  not  agree  with  Stratton's 
assumption  for  the  ion  gyro  radius.  Nevertheless,  the  value  of  that 
Stratton  obtained  from  these  assumptions  are  within  the  same  order  of 
magnitude  as  the  values  obtained  from  the  computer  simulation. 


I.  Discussion 

From  the  collisionless  model  of  the  direct  electromagnetic  ac¬ 
celerator,  one  can  deduce  (l)  that  the  energy  of  the  accelerated  ions, 
in  electron  volts,  cannot  exceed  the  voltage  applied  to  the  accelerator, 
and  (2)  that  the  optimum  performance  would  occur  with  no  applied  mag¬ 
netic  field.  These  are  in  contrast  to  the  electric  accelerators  op¬ 
erating  in  the  MPD  regime  where  experimental  measurements  have  shown 
that  the  accelerated  ions  can  receive  an  energy  in  electron  volts 
greater  than  the  voltage  applied  to  the  accelerator  and  wherein  an 
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applied  magnetic  field  can  improve  the  performance  .  When  charged 
particles  receive  energy  only  from  steady  fields,  they  cannot  obtain  a 
higher  energy  than  that  of  falling  through  the  maximum  difference  in 
potential.  For  the  colliaionless  model,  the  energy  of  the  ions  leaving 
the  accelerator  will  equal  in  electron  volts  the  difference  in  potential 
between  the  point  of  ionization  and  the  potential  outside  of  the  ac¬ 
celerator  exhaust.  The  maximum  thrust  will  occur  when  all  of  the  energy 
the  ion  receives  falling  through  the  potential  is  converted  to  axial 
velocity.  The  production  of  an  azimuthal  velocity  would  merely  shift 
part  of  the  axially  directed  energy  to  a  rotational  energy.  Since  a 
steady  axial  or  radial  applied  magnetic  field  produces  rotation  about 
the  axis  but  cannot  add  energy  to  the  ions,  the  rotational  energy  would 
be  subtracted  from  the  axially  directed  energy  resulting  in  a  reduced 
thrust . 

The  coulomb  cross  section  for  high  energy  collisions  of  charged 
particles  with  ions  becomes  less  than  the  hard  sphere  cross  section 

of  the  ions  if  the  ions  have  electron  shells.  For  example,  the  ion- 

4 

ion  cross  section  for  hydrogen  ions  with  10  ev  energy  is  of  the  order 
-20  2 

of  10  cm  whereas  the  hard  sphere  cross  section  for  ions  with  electron 

-!«=  p 

shells  is  of  the  order  of  10  ''  cm  '.  Therefore,  although  the  hydrogen 

ion  mean  free  paths  may  be  largo  compared  to  the  size  of  the  accelerator 
for  the  high  specific  impulse  accelerators,  the  mean  free  path  of 

other  ions  would  be  of  the  order  of  the  size  of  the  accelerators  at 

lc  3 

ion  number  densities  of  10  'em  ,  Therefore,  the  collisionless  model 
would  not  apply  rigorously  for  gases  other  than  hydroger.  Nevertheless, 
the  ions  would  still  receive  most  of  their  energy  directly  from  the 
electric  field,  and  the  collisionlesr;  model  may  still  predict  the 
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accelerator  characteristics  satisfactorily.  (The  influence  of  the 
atomic  weight  and  ion  charge  on  the  characteristics  in  the  collisionless 
model  will  be  discusser,  in  the  next  chapter  entitled  "Performance 
Predictions  for  the  Direct  Electromagnetic  Accelerator") 

The  voltage-current  characteristics  given  by  the  computer  sim¬ 
ulations  are  shown  in  figure  15.  It  was  assumed  that  large  axial  electric 
fields  at  the  cathode  would  induce  very’  large  currents  from  the  cathode, 
simulating  high  electrical  conductivity  in  and  near  the  cathode.  The 
simulation  was  started  with  an  arbitrarily  low  current  and,  as  the  sim¬ 
ulation  progressed  in  time,  the  current  increased  to  an  asymptotic  value. 
As  the  current  increased,  the  self-pinch  caused  the  axial  electric 
field  near  the  cathode  to  decrease  to  a  small  value  as  shown  in  figure 
12  so  that  the  asymptotic  currents  do  not  become  excessively  large. 

Since  there  are  no  particle-particle  collisions  in  the  simulation,  the 
electrical  conductivity  is  essentially  infinite.  The  radial  electron 

current  throughout  the  accelerator,  however,  is  limited  by  the  fact  that 

— ♦  —♦ 

the  electric  fields  are  balanced  by  the  v  X  B  force  on  the  electrons. 

e 

The  ion  currents  are  limited  by  inertia  forces  in  addition  to  the  partial 
balancing  of  the  electric  field  with  the  v^  X  B  force. 

Although  the  computer  simulation  was  developed  primarily  for 
comparison  with  the  qualitative  assumptions  made  for  the  simplified 
models,  it  may  be  used  to  simulate  the  plasma  flows  in  accelerators 
of  various  types  to  determine  which  configurations  will  be  the  most 
efficient.  For  example,  introducing  the  gas  at  the  anode  and  forcing 
it  to  be  ionized  in  the  anode  sheath  would  force  the  ions  to  fall 
through  a  greater  potential  than  if  the  gas  were  ionized  nearer  to 
the  cathode.  Therefore,  an  accelerator  with  the  gas  introduced  at  the 
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anode  should  produce  higher  thrust  than  one  with  the  gas  introduced 
near  the  cathode.  A  simulation  was  performed  with  the  ionization 
region  assumed  to  be  near  the  anode  rather  than  distributed  according 
to  figure  9>  and  an  increase  in  thrust  was  noted. 

The  computer  simulation  could  be  extended  to  include  cases  with 
applied  magnetic  fields.  This  extension  would  require  additional 
computer  storage  space  in  order  to  store  the  values  of  the  radial  and 
axial  components  of  the  magnetic  field  and  to  store  the  azimuthal 
velocities  of  the  particles.  The  azimuthal  force  equation  would  have 
to  be  added  in  order  to  calculate  these  azimuthal  velocities,  and  the 
azimuthal  currents  would  have  to  be  computed  from  these  velocities  in 
order  to  calculate  the  induced  axial  and  radial  magnetic  fields.  These 
additional  calculations  would  increase  the  computation  time  by  an  es¬ 
timated  20^. 

The  values  of  velocity,  thrust,  and  voltage  obtained  with  the 
computer  simulation  of  Cheng's  accelerator  agree  with  the  experimental 
measurements  of  Cheng  to  within  the  accuracy  of  the  experimental  measure¬ 
ments. 

The  current  distribution  within  the  simulation  is  similar  to  the 

4 

current  distributions  measured  by  Clark  and  Jahn  in  that  current  locps 
protrude  downstream  of  the  anode.  However,  further  comparison  is  not 
yet  possible  because  Clark  and  Jahn  have  not  yet  reported  values  of 
velocity,  thrust,  or  voltage  within  their  high  current  accelerator 
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V.  PERFORMANCE  PREDICTIONS  OF  THE  DIRECT  ELECTROMAGNETIC  ACCELERATOR 

A.  Dimensionless  Parameters 

When  the  equations  used  in  the  computer  simulation  program  (given 
in  the  appendix)  are  made  dimensionless,  the  significant  parameters  that 
emerge  are 


where 

Z  =  the  number  of  charges  on  each  ion 
e  =  the  charge  of  a  singly  ionized  ion,  coulombs 
=  mass  of  the  ion,  kg 

$  =  potential  applied  across  the  electrodes,  V 

At  =  incremental  time  step  between  calculations,  s 

L  =  characteristic  dimension  of  the  accelerator,  m 

2 

H  =  magnetic  permeability,  N/A 
I  =  current ,  A 
m  =  mass  flow  rate,  kg/s 

The  parameters  above  are  fixed  for  each  numerical  simulation,  but  the 
simulation  applies  to  any  combination  of  applied  voltages,  currents, 
mass  flow  rates,  accelerator  sizes,  and  time  scales  that  give  the  same 
dimensionless  parameters.  For  example,  one  simulation  that  represents 
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an  applied  potential  of  1000  volts,  a  current  of  20,000  amperes,  and  a 
hydrogen  flow  rate  of  .2  g/s  also  simulates  the  same  accelerator  with 
an  applied  potential  of  9000  volts,  a  current  of  60,000  amperes,  and 
a  hydrogen  flow  rate  of  .6  g/s.  The  movie  of  the  simulation  also 
represents  either  case,  but  the  incremental  time  step  for  the  second 
case  Is  1/3  the  incremental  time  step  for  the  first  case  solhe  time 
dilitation  of  the  movie  is  different  for  the  two  cases. 

From  these  dimensionless  parameters,  the  variation  of  current, 
mass  flow  rate,  specific  impulse,  and  thrust  as  a  function  of  the 
applied  voltage,  the  atomic  weight,  and  the  number  of  charges  of  each 
ion  for  a  given  accelerator  are  as  follows: 

I  «  (tA/Z)1/2 

i  «  »1/Vz)3/a 

1/2 

specific  impulse  «  ($Z/A) 

thrust  =  $A/Z 

where  A  =  atomic  weight 

/ 

B.  Predicted  Performance  Charts 

The  performance  obtained  from  the  computer  simulation  for  Clark 
and  Jahn's  accelerator  and  for  Cheng's  accelerator  are  shown  in  figures 
16  and  17.  For  voltages  below  the  limit  marked  "collisional  plasma" 
on  figure  1 6,  the  collisionless  approximation  is  not  expected  to  be 
realistic,  and  for  voltages  above  the  limit  marked  "relativistic  electrons" 
one  would  have  to  include  the  relativistic  effects  in  the  computer  sim¬ 


ulation  . 


The  values  given  on  figure  l6  are  independent  of  the  size,  but 
the  thrust  density  shown  on  figure  17  depends  on  the  cross  sectional 

area  of  the  accelerator.  The  cross  sectional  area  of  Cheng's  ac- 

2 

celerator  is  approximately  10  cm  whereas  Clark  and  Jahn's  is  approx- 

2 

imately  100  cm  .  This  difference  in  cross  sectional  area  is  the  primary 
reason  the  curve  for  Clark  and  Jthn's  accelerator  is  so  far  below  the 
curve  for  Cheng's  accelerator. 

When  the  assumed  location  for  the  region  of  ionization  was  varied 
with  the  applied  voltage  held  constant,  the  current,  mass  flow  rate, 
velocity,  thrust,  and  specific  impulse  varied  less  than  a  factor  of 
two.  Therefore,  these  predictions  are  considered  to  be  good  enough 
for  order  of  magnitude  estimates  despite  the  fact  that  the  location 
for  the  region  of  ionization  is  unknown  and  had  to  be  assumed. 


VI .  CONCLUSIONS 


The  ions  within  the  very  high  specific  impulse  accelerators  (such 
as  Cheng's  coaxial  accelerator)  are  accelerated  directly  by  the  electro¬ 
magnetic  fields  rather  than  by  collisions  with  the  electrons.  For  this 
regime  of  direct  electromagnetic  acceleration,  the  present.  MPD  arcjet 
theories  are  inapplicable. 

A  collisionless  plasma  model  is  recommended  to  represent  the  region 
where  the  ions  are  receiving  their  energy  directly  from  the  electro¬ 
magnetic  fields.  The  Lagrangian  computer  simulation  technique  appears 
to  simulate  satisfactorily  the  plasma  within  the  collisionless  region 
of  these  direct  electromagnetic  accelerators.  The  values  of  velocity, 
thrust,  and  voltage  obtained  with  the  computer  simulation  of  Cheng’s 
accelerator  agree  with  the  experimental  measurements  of  Cheng  to  within 
the  accuracy  of  tl»?  experimental  measurements. 

Results  from  the  computer  simulation  show  that  the  simple  par¬ 
ticle  models  proposed  by  Stratton  and  the  author  contain  the  essential 
features  for  describing  the  acceleration  mechanism  within  direct  electro¬ 
magnetic  accelerators.  Order  of  magnitude  estimates  of  the  performance 
characteristics  of  the  direct  electromagnetic  accelerator  obtained  from 
the  computer  simulation  are  in  agreement  with  the  predictions  of 
Stratton’s  collisionless  particle  model  despite  the  fact  that  some  of 
Stratton's  assumptions  are  not  corroborated  by  the  computer  simulation. 

Results  of  the  computer  simulation  also  illustrate  that  the 
p^E  force  term  does  not  strictly  counteract  the  p^E  force  term  even 
though  the  plasma  is  nearly  neutral.  Rather,  the  pgE  force  is  balanced 
by  the  xB  force,  and  the  p^E  force  is  the  dominant  force  for  adding 
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energy  to  the  ions.  In  addition,  even  though  the  electrical  con¬ 
ductivity  is  essentially  infinite  because  no  particle-particle  en¬ 
counters  are  considered,  the  plasma  is  not  frozen  to  the  magnetic 
field  lines  because  of  the  large  ion  inertia  forces. 

The  direct  electromagnetic  accelerator  should  operate  most 
efficiently  with  no  applied  magnetic  field  and  with  a  mass  flow 
proportional  to  the  current.  The  specific  impulse  should  increase  in 
proportion  to  the  current,  and  the  voltage  and  thrust  should  increase 
with  the  square  of  th9  current. 
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APPENDIX.  COMPUTER  PROGRAMS  FOR  SIMULATING  THE 
DIRECT  ELECTROMAGNETIC  ACCELERATORS 


These  computer  programs  were  written  in  FORTRAN  IV  and  were 
run  on  the  IBM  7094  DCS  system  at  Ames  Research  Center,  NASA.  The 
programs  are  located  as  follows: 

PROGRAM  PAGE 

1.  Program  For  the  Simulation  of  Clark  and  Jahn's 
Accelerator 

a.  Main  program . 51 

/ 

/ 

b.  Subroutine  INITL.  ..  ." . 55 

c.  Subroutine  INITL  (for  a  continuation  run)  .  .  .  60 

d.  Subroutine  FIELDS  .  65 

e.  Subroutine  POISS . 73 

f.  Subroutine  TRAJ . 75 

g.  Subroutine  SOURCE  .  . 84 

h.  Subroutine  OUTPT . 90 

1.  Subroutine  OTTPT . 93 
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PROGRAM 


PAGE 


2.  Program  for  the  Simulation  of  Cheng's  Accelerator 


a.  Main  program . • . 97 

b.  Subroutine  INITL . 101 

c.  Subroutine  INITL  (for  a  continuation  run)  .  .  .  106 

d.  Subroutine  FIELDS  .  Ill 

e.  Subroutine  POISS . 73 

f .  Subroutine  MAT . 119 

g.  Subroutine  TRAJ  . . 121 

h.  Subroutine  SOURCE . 130 

1.  Subroutine  OUTPT.  . . 136 

j.  Subroutine  OTTPT . 139 
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180  CONTINUE 
LAST  =  1 
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* EZ (3*15) •  ERi 5* 15  ) 

REAL  T»OT.MASSR»NAMPL.EPSAMP,REFLC» VOLTS* EX #DR*DZ*PI *Q*MUZESO*EPSO 
•  MASSI . MASSE. RDS » LENG . AMP5E , AMPS  I »QNET 
real  EXTRA1*EXTRA2* EXTRAS 


INTEGER  NsE<5,15>,  NS*(5,15)»  NDELl<400>*  NDELE(4qo> 
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DO  100  N*1 *NI T 
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***  ELECTRON  TRA J~CTOR I ES  *** 
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***  BOUNDARY  CONDITIONS  FOR  THE  PARTICLES  *** 
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CHANGING  back  to  arrays  *** 
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Figure  1.  Performance  of  Various  Propulsion  Devices  as  Depicted  by 
Clark  and  Jahn1 
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Figure  2#  A  Conventional  Constricted  Arc jet 
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Figure  4.  Cross  Section  of  Clark  and  Jahn's  Pulsed  Accelerator 


Figure  7„  Mean  Free  Paths  of  Hydrogen  Ions  at  Ion  Number  Densities  of 
10  ions/cc 
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Figure  8.  A  Typical  MPD  Accelerator  (From  Jahn^) 
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CHENG'S  ACCELERATOR  CLARK  a  JAHN'S  ACCELERATOR 


ASSUMED  REGION 


Figure  9.  Accelerator  Models  Used  in  the  Computer  Simulations 
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CATHODE-1  INSULATOR- 


Figure  10.  Mesh  Configuration  Within  the  Models  of  the  Accelerators 
Used  for  the  Computer  Simulation 
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LEGEND 

The  representative  particles  simulating  the  plasma;  the  large  dots 
represent  ions  and  the  small  dots  represent  electrons. 

The  reference  plane  of  constant  aximuthal  position  showing  the 
radial  and  axial  positions  within  the  accelerator. 

The  same  reference  plane  translated  to  a  new  position  in  order 
to  plot  the  electric  potential  and  the  magnetic  field  as  a 
function  of  radial  and  axial  position.  For  each  point  in  the 
reference  plane,  the  vertical  height  to  the  surface  plotted 
represents  the  magnitude  of  the  electric  potential  or  of  the 
magnetic  field  for  that  radial  and  axial  position. 

The  number  of  representative  electrons  and  ions  currently  being 
used  to  simulate  the  plasma. 

A  clock  measuring  the  time  in  units  of  transit  times.  A  transit 
time  is  the  time  that  would  be  required  for  an  ion  to  travel 
the  length  of  the  accelerator  if  it  had  an  energy  in  electron 
volts  equal  to  the  potential  difference  between  the  electrodes  in 
volts . 

F.)  The  number  of  movie  frames  from  the  start  of  each  case. 


Figure  11.  Sample  Frame  of  the  Movie  of  the  Plasma  Simulation 
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Figure  12.  Modification  of  the  Applied  Potentials  by  the  Plasma 


Figure 


Typical  Particle  Trajectories  Observed  in  the  Computer 
Simulation. 


AXIAL  POSITION 


Figure  l4.  Electron  and  Ion  Trajectories  In  a  Uniform  Electric  Field  And 
A  Magnetic  Field  That  Decays  in  the  Axial  Direction 


CURRENT,  A 


Figure  15.  Voltage-Current  Characteristic.  Obtained  from  the  Computer 
Simulation 
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Figure  16.  Performance  Predictions  Obtained  from  the  Computer  Slamlations 


cu 

E 


>- 

t 

c/> 

z 

Ui 

o 

}“ 

</> 

3 

a: 


I07  r- 

I06  - 
I05  - 
I04  - 
I03  - 
I02  - 
10  - 
I  ^ 


- , - r 

COMPUTER  SIMULATION  OF- 

cheng^s  accelerator 


CHENG'S  EXPERIMENTAL 
DATA - v 


CNVEITIMAl 

MCJET 


-COMPUTER  SIMULATION  OF  _| 
CLARK  ft  JAHN'S  ACCELERATOR 


,/////, 

PULSED 


I0‘ 


10' 


10' 


I06 


SPECIFIC  IMPULSE,  sec 


Figure  17.  Performance  Predictions  Obtained  from  the  Computer  Simulations 
Compared  with  Cheng's  Experimental  Data  and  With  the  Perfor¬ 
mance  of  Various  Propulsion  Devices  as  Depicted  by  Clark  and 
Jahn* 


